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It is usual to study confinement via quantum chromodynamics (QCD) alone. The deconfinement 
transition of the pure gauge theory (i.e. with static quarks) is then characterized by the breaking 
of center symmetry. Center vortices offer an intuitive and quantitative description of the tran- 
sition. Dynamical quarks explicitly break center symmetry, and the phase transition becomes a 
crossover. However, it may be misleading to study QCD in isolation. Quarks also carry fractional 
electric charge. This bestows the Standard Model with a global center symmetry that combines 
color center phases with an appropriate electromagnetic phase. Is this symmetry relevant to con- 
finement? We begin our investigation by studying a 2-coIor model of QCD with half-integer 
electrically charged quarks. 
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1. Introduction 

The finite temperature deconfinement transition of pure SU(A'^) gauge theory is intimateiy tied 
to center symmetry. The expectation vaiue of the Polyakov loop {P) jumps from zero to a finite 
value at Tc, thereby picking a center sector and signaling deconfinement along with the sponta- 
neous breaking of center symmetry. The transition is well described by the dynamics of spacelike 
center vortices, whose proliferation at low temperatures disorders the Polyakov loop and leads to 
confinement [1]. They carry center flux and separate regions where the Polyakov loop differs by a 
phase z € Zj^. The suppression of center vortices at high temperature coincides with the ordering 
of the Polyakov loop, and their free energy offers an elegant dual order parameter for the transition 
[2, 3]. Deconfinement mirrors the order-disorder transition of a simple Zat spin model [4], with 
vortices playing the role of spin interfaces. When both phase transitions are of second order their 
behavior is universal. We have fruitfully exploited this for SU(2) and SU(3) in 2-1-1 dimensions, 
where many exact results are available from the corresponding spin models [5,6]. 

Center symmetry is explicitly broken by the introduction of quarks, however, and this picture 
is lost. Their ordering effect is easily seen via a hopping expansion of the fermion determinant. 
One obtains Polyakov loop terms that favour the {P) = 1 center sector in the same way that an 
external magnetic field aligns spins in a spin system. 

Dynamical quarks also make the role of center vortices unclear. Vortex sheets, which were 
freely moved in the pure gauge ensemble, become observable. Quarks pick up the center flux as a 
phase when they encircle a vortex. But what about the quarks' fractional electric charge? Phases 
from the center of the color gauge group may be compensated by an appropriate electromagnetic 
contribution. Phase cancellation renders a combined center vortex and electromagnetic Dirac string 
unobservable. This reflects a hidden global gauge symmetry of the Standard Model. It combines 
the centers of each group factor and is therefore an extension of center symmetry. 

We would like to study the relevance of this symmetry to confinement and the phase structure 
of the Standard Model. As a first step we consider 2-color QCD with the addition of quarks carrying 
fractional electric charge. Here we concentrate on the possibility that the quarks' electromagnetic 
coupling leads to disorder in the color fields, akin to placing a spin system in a random external 
field. 

2. The hidden symmetry 

When subjected to a nontrivial SU(3) center gauge transformation, or upon encircling a center 
vortex, quarks pick up a phase g^i^'^^/^ ^ 23. Given their electric charges, 2 = |e or —\e, an 
additional electromagnetic factor g^'^'^'^S/e enough to exactly cancel the color contribution. This 
is in fact the phase that a particle with electric charge Q receives when it encircles the Dirac string 
of a monopole with minimal magnetic charge with respect to e. Colorless particles with integer 
electric charge are blind to both phases by default. Therefore, the pair of combined transformations 

^^anl^^^anQ/e^^ ^^-ilnl^ ^^-ilnQj ^ 3^(3) X U(l)„„ (2.1) 

act trivially on all particles in the Standard Model. Together with the identity transformation they 
form a global Z3 symmetry. The symmetry doubles when we include the weak interactions. Electric 
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charge Q is related to hypercharge Y and the third component of weak isospin tT,hy Q/ e = tj, + Y /2. 
Since e'^^'^ represents —1 G SU(2), the electromagnetic phase g'^'^S/e given by {—\,e'^^) G 
SU(2)xU(l)}'. The discrete symmetry group is therefore generated by 

{e'^''/\-l,e"'^) G SU(3) x SU(2) x U(l)}.. (2.2) 

This gives six elements, so the 'true' or minimal gauge group of the Standard Model is 
SU(3)xSU(2)xU(1)f/Z6. See [7] for a review. 

This Zg symmetry is very important in the context of grand unified theories. Without it, the 
Standard Model gauge group would be too large to fit in SU(5) C SO(IO) . Moreover, it endows 
the model with non-trivial topology that allows for defects carrying both color and electromagnetic 
flux. These are, namely, combined center vortex and electromagnetic Dirac strings/surfaces. Given 
the importance of center vortices in pure SU(A'^) gauge theory, it is natural to ask if the discrete 
symmetry and combined vortices are relevant to the confinement of dynamical quarks.^ This type 
of defect unification has already been shown to have a large impact on the phase diagram of a toy 
U(l)xU(l) model [13, 14]. 

3. Toy model 

Lattice simulations of the full Standard Model gauge group with dynamical quarks are beyond 
our reach. As such, we start with a toy model of 2-color QCD plus electromagnetism with 2 flavors 
of Wilson fermions in 3+1 dimensions. Our 'up' and 'down' quarks have fractional chai^ge zt^e 
relative to the U(l)e,„ gauge action. The weak interactions are discarded. Nevertheless, we have 
a chance to see what interesting physics may come from the quarks' non-trivial coupling to both 
photons and gluons. The lattice action is 

S=-Y. (^ReTr[/p + j8,^cos0p)+S/,w, (3.1) 

plaq. V ^ / 

where Sf^w is the usual Wilson fermion action with the distinction that the parallel transporters for 
quarks are products of an SU(2) color matrix and a U(l)em phase. They are of the form, 

[/^(x)e'^''W/2, Uf,{x) G SU(2), 0^(x) G (-271,271]. (3.2) 

The SU(2) plaquettes Up and U(l) plaquette angles dp are formed from and 6^ in the usual way. 
The fractional charge of the quarks is encoded in the fact that their parallel transporters contain half 
the U(l) angle relative to the 6^ that appear in the plaquette angle dp. For instance, what the quarks 
see as an e'^f /2 = _i electromagnetic link appears as a e'^^' = +1 link in the gauge action. Since 
the 'size' of the compact U(l) is determined by the quarks' electric charge, the range of is 
chosen such that we integrate over all possible electromagnetic transporters for the quarks. 

The gauge fields are blind to global SU(2) center and U(l) transformations, so we would have 
a global Z2xU(l) center symmetry in the absence of fundamental fields. It is clear from Eq. (3.2) 
that a combined color center phase — 1 G SU(2) and electromagnetic phase of — 1 acts identically 

^The Zg symmetry was recently studied on the lattice by Bakker, Veselov and Zubkov [8]. There they were more 
interested, however, in consequences for electroweak and Higgs physics. 
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on the quarks. Thus a Z2 subgroup of this symmetry survives their introduction, and the minimal 
gauge symmetry of our model is SU(2)xU(l)em/22- ^2 plays the role of the hidden symmetry of 
the Standard Model. 

The loops in a hopping expansion of the fermion determinant are now constructed from both 
color and electromagnetic parallel transporters. For Nt = 4 time slices, the leading contribution 
comes from plaquette and Polyakov loop terms. Arranged in an effective action, they give 

5/,eff = -16fc4 ( £cos^-ReTr[/p + 8X;cosPe/2-ReP,„, j + ..., (3.3) 

\plaq. X / 

where k is the hopping parameter. Note that the U(l) contribution to the plaquette term is formed 
from U(l) link angles as seen by the quarks, i.e. e'^f/^. P^oi denotes the traced SU(2) Polyakov 
loop, while Pq/2 denotes the U(l) Polyakov loop angle, Pe/2ix) = L|1o' ^t{t,x)/2. 

In the absence of electromagnetism, the Polyakov loop term has the same effect as coupling 
an external ordering magnetic field to a spin system. It picks an orientation. The U(l) coupling 
alters this picture, however, since disorder in the quarks' U(l) links promotes disorder in the SU(2) 
hnks. If cosPe/2 is disordered, then the Polyakov loop term from the hopping expansion will be 
minimized only if P^.gi is also disordered. The effect on P^.gi may be compared to placing a spin 
system in a random magnetic field. For instance, consider an Ising model with the Hamiltonian 

Jif = -j'£siSj-h'£hiSi, (3.4) 

('•.;) 

where 5, = ±1 are spins and the direction of the magnetic field /?,• = ±1 is site dependent and 
random. Such models have very interesting phase structures. Order-disorder phase transitions are 
still possible, but their existence and type depend on the magnitude of the applied random field [9]. 

More realistically, imagine an Ising model in a fluctuating external field. If the fluctuations 
are large they will hinder the ordering of spins, thus preserving the global Z2 symmetry. This is 
reminiscent of the Peccei-Quinn mechanism in QCD [10], which replaces the 6 parameter of the 
CP violating term by an axion field. Dynamical preservation of CP symmetry is then achieved if 
the fluctuations of the axion field average to zero. 

In the case of spins, we can model the interaction with a dynamical field by including addi- 
tional U(l) variables at each site. This gives a so-called XY-Ising model [11]. Indeed, the typical 
coupling of U(l) and Ising energy densities in these systems is similar to the mixed plaquette term 
from Eq. (3.3). As for spins in a static random magnetic field, the coupling in an XY-Ising model 
leads to very novel critical behavior. We may expect the same for our gauge theory. 

Compact QED has a disordered, confining phase for jSg,,, < 1. Here the U(l) should have a 
large effect. However, this is not a physically relevant regime since it confines all electric charges. 
We would like to confine quarks but not integer charged particles such as electrons. The interesting 
question is whether the U(l) has an influence on the color dynamics in the Coulomb phase for 
integer charges, in which our U(l) links are ordered with respect to the pure gauge action. Now 
the fractional electric chai^ge of quarks becomes important. The phases 6^ and 6^ + 27i are distin- 
guished by the quarks, but not by the gauge action. As we cross fiem ~ 1, the U(l) links should 
become ordered for integer charges, but there is still room for Z2 disorder as seen by quarks. Does 
this lead to disorder in the color links? 
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Figure 1: Volume averaged SU(2) Polyakov loop in the presence of quarks with fractional electric charge. 
On an 8^^ x 4 lattice with K ~ 0.15. The 'unquenched' and 'pure' results are without electromagnetism. The 
left plot is from a hot start, the right from a cold start. Results with fc = 0.175 (not shown) are similar. 
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3.1 Dynamical simulations 

We checked this by measuring the magnitudes of volume averaged Polyakov loops in our 
SU(2)xU(l)/Z2 toy model with dynamical quarks. We update via Hybrid Monte Carlo.^ For 
initital testing, we restricted ourselves to 8^ x 4 lattices. We have yet to calculate the mass scale 
and critical in our model, so we use the 2-color results of Skullerud et al. for guidance [12]. At 
zero temperature they find a lattice spacing of c?^' w 690 MeV for jS^o/ = 1.8, K = 0.178. From 
this we chose K = 0.15, 0.175 for 'heavy' and 'hght' simulations, respectively. In Fig. 1 we plot 
the results for the SU(2) Polyakov loop measured near the pure SU(2) phase transition. Various 
U(l) coupUngs were used, from the U(l) confined phase (j3e„, = 0.8) to deep into the Coulomb 
phase for integer charges (j8e,„ = 8). We verified that the U(l) Polyakov loop for integer charges 
(i.e. constructed from e'^' links) jumps from zero to one in the vicinity of Pem = l^ which signals 
that the compact QED phase transition remains essentially unchanged. 

Using a random, hot start for both the color and U(l) fields, the pure gauge result is reproduced 
for the SU(2) Polyakov loop, despite the fact that our model contains dynamical quarks (see Fig. 1 
left). This is the case for both our 'heavy' and 'light' simulations. We find the same result using 
a cold start with respect to the U(l) gauge action but with Z2 disorder in the U(l) links as seen by 
quarks, i.e. = I, e'^f = ±1. The U(l) gauge action is unable to dynamically remove the Z2 
disorder and it becomes frozen in. This disorders the SU(2) Unks via the coupling to quarks and 
one recovers the result for the pure gauge theory. 

Things are more complicated when we simulate from a completely ordered start (U^ = 1, 
_ jj^g quenched SU(2) behavior persists into the Coulomb phase for integer electric 
charges, but not indefinitely. We observe a transition to the unquenched SU(2) result near jS^m = 1.3 

^From an algorithmic point of view we use 2 degenerate (i.e. equal charge) species of quarks. However, owing to 
the pseudoreality of the fundamental representation of SU(2), this is equivalent to a system with i^e quarks. You can 
convince yourself of this by playing with the charge conjugation operator and color basis. 
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Figure 2: Volume average SU(2) Polyakov loop, where the effect of electrically charged quarks is included 
up to 0'{k'^) in the hopping expansion. On an 24^ x 4 lattice with K = 0.175. The 'unquenched' results are 
without electromagnetism. The left plot is from a hot start, the right from a cold start. 



(see Fig. 1 right). Fluctuations of the link angles are suppressed so strongly that the quarks are 
unable to generate Z2 disorder. In this case, the order from our cold start is frozen in and effectively 
breaks the Z2 symmetry. 

3.2 Hopping expansion simulations 

Before performing full dynamical simulations on larger lattices, we would like to see how far 
our intuition from the hopping expansion can take us. As such we have run simulations for Nt = 4 
where the effect of quarks are included up to ^{k^). This means that we include the terms in Eq. 
(3.3) in addition to the pure gauge actions for SU(2) and U(l). We simulate using K = 0.175 on 
24^ X 4 lattices, once more from both hot and cold starts. The results are shown in Fig. 2. As for 
the fully dynamical case, Z2 disorder from a hot start leads to quenched behavior for the SU(2) 
Polyakov loop. The red points are from simulations of our hopping expansion effective action for 
SU(2) alone. From a cold start we observe a transition towards this curve near j8e,„ > 1 .3. We have 
crosschecked this behavior using both HMC and a mixture of SU(2) heatbath plus U(l) Metropolis 
updates. This suggests that the retention of Z2 disorder from a hot start is a feature of the model 
and is not an artifact of our HMC algorithm for fully dynamical quarks. 

Since the hopping expansion captures the essential features of the model, it can be exploited 
to better understand how the fractional electric charge of quarks influences the color dynamics. 
For instance, the Polyakov loop term has the dominant effect for unquenched SU(2) without elec- 
tromagnetism. The plaquette term merely shifts Pcoi- With the coupling to U(l), however, the 
plaquette contribution from the hopping expansion becomes relevant. It also conveys disorder 
between the color and electromagnetic links. By simulating with and without each term we can 
determine their relative importance. This may help us to construct an effective spin model. Via the 
hopping expansion we can also experiment with large values of K. From a hot start one may then 
expect the color and U(l) links to align with respect to each other but not with respect to the gauge 
action. These points will be addressed in the future. 
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4. Discussion 

The results for our SU(2)xU(l)/Z2 model demonstrate that the fractional electric charge of 
quarks relative to a compact U(l)em gauge action may have a non-perturbative effect on the color 
sector. It is related to a discrete global symmetry that combines electromagnetic and color phase 
factors. Although the electromagnetic phases are not distinguished by the U(l) gauge action, dis- 
order amongst them leads back to the color degrees of freedom via the coupling to quarks. This 
is enough to recover quenched behavior for the SU(2) Polyakov loop, even far into the Coulomb 
phase for integer charges. 

Since the Standard Model contains such a global symmetry, it may be misleading to study the 
group factors SU(3)co; and \]{\)em separately. It is tempting to speculate about how the fractional 
charge of quarks may alter the 'Colombia plot' picture of confinement [15]. For QCD with zero 
chemical potential, a first order transition persists for heavy quarks. If the nontrivial coupling to 
U(l)e„, is able to disorder the color fields, this may extend the first order region to lighter quarks 
and perhaps sharpen the crossover near the physical quark masses. Note, however, that our results 
depend on the initial conditions of the simulations. In order for the U(l)e„, to have an influence, 
we need a scenario in which disorder with respect to the electromagnetic phases seen by quai^ks is 
frozen in when we reach the perturbative regime for integer electric charges. This leads naturally 
to a consideration of grand unified theories. It also prompts us to consider the type of field con- 
figurations that provide the disorder. We expect these to be vortices that carry both color and elec- 
tromagnetic flux. Generalizing 't Hooft's twisted boundary conditions [16] to our SU(2)xU(l)/Z2 
model, we will be able to study combined vortices in prescence of dynamical quarks. 
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